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Abstract 


The ion motions in a Penning trap have been studied in detail in the presence of azimuthal dipolar and quadrupolar radio- 
frequency excitations and buffer gas cooling. The numerical solutions by using the Runge-Kutta method and thus the pictures of 
the ion trajectories in the trap have been obtained for different cases and summarized in graphical form. For the recentering of the 
ion of interest and to perform the purification of the ion species, one has to set a reasonable buffer gas pressure in the trap and apply 


azimuthal quadrupolar excitation at frequency w,f = We. 
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1. Introduction 


Penning traps have become very accurate tools for mass de- 
termination both on stable and unstable isotopes. There are 
many Penning traps which are in operation or under construc- 
tion all over the world, such as ISOLTRAP [1] in CERN, 
SHIPTRAP [2] in GSI, LEBIT [3] in MSU, CPT [4] in ANL, 
JYFLTRAP in JYFL [5] and so on. At the Institute of Mod- 
ern Physics, Chinese Academy of Sciences, the LPT (Lanzhou 
Penning Trap) is also under construction [6]. 

In an ideal Penning Trap, an ion is confined by the combi- 
nation of an electrostatic quadruploar field and a homogeneous 
magnetic fields, and its motion is a superposition of three eigen- 
motions, an axial oscillation (z) with frequency w; and two ra- 
dial motions commonly referred as reduced cyclotron (+) and 
magnetron motions (—) with frequency w+ and w_, respectively. 
To measure the ion’s mass with high precision, one of the meth- 
ods is to perform buffer gas cooling and radio-frequency (rf) ex- 
citations to remove unwanted ions from nuclear reactions and 
other sources. This method has been used in many Penning 
traps. 

The ion motion can be driven by oscillating electric fields 
and this in general results in a change of the amplitudes of the 
eigenmotions. The effect of the driving field on the ion mo- 
tion depends on the multipolarity of the field and its frequency. 
A dipole field at one of the eigenfrequencies can be used to 
increase the amplitude of the corresponding eigenmotion, and 
this dipolar excitation at w,f = w- is generally a significant step 
to remove the unwanted ions. If ion is excited by a quadrupo- 
lar field with a frequency w,f = we, where we is the cyclotron 
frequency of the ion, the magnetron motion and cyclotron mo- 
tion will be continuously converted into each other, and this is 
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always used to determine the true value of we. If all three eigen- 
frequencies are measured, the invariance theorem of Brown and 
Gabrielse [7] can be used. 


Cooling of stored ions results in a reduction of the motional 
amplitudes, and the cooled ions can be trapped in a much 
smaller volume and thus probe less of the imperfections in the 
trapping electric and magnetic fields. The technique of buffer 
gas cooling is commonly applied for radioactive ions stored in 
a Penning trap, and noble gases are typically used because of 
their high ionization potential and thus minimum losses due to 
charge exchange. 


So an ion confined in a Penning trap will experience the 
forces from both the electrostatic field and the magnetic field, 
and the effects of rf excitation and buffer gas. The equations of 
ion motion become very complicated and it is very difficult to 
solve them analytically without any approximation. 


Starting from analytical solution, Bollen et al. [8] described 
ion trajectories under different excitation scenarios qualitatively 
and identified the important causes of uncertainty in high- 
precision mass measurements of heavy ions. Savard et al. [9] 
presented ion trajectories in Penning traps in the presence of 
buffer gas from Runge-Kutta integration of the relevant equa- 
tions of motion and demonstrated experimentally the effect of 
buffer gas cooling. König et al. [10] investigated the ion motion 
in the presence of an azimuthal quadrupolar rf field and buffer 
gas cooling and excellent agreement was observed between the- 
oretical results and experimental data. 


In this paper we solve those equations of ion motion numer- 
ically by using Runge-Kutta method, obtain the ion trajectories 
in the trap for many different cases and summarize the ion tra- 
jectories in graphical form to help us to understand the physical 
picture in a Penning trap. 
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Figure 1: Schematic layout of a typical Penning trap 


2. Dynamical equations of ion motion 


Figure | shows the schematic layout of a typical Penning 
trap. In an ideal Penning trap a charged particle with a mass 
of m and a charge of q is confined by the combination of a 
homogeneous magnetic field B= Bê, and an axial quadrupo- 
lar potential ®(z, p) = oo (22 - 5p”), where Uo is the applied 
trapping voltage between the ring electrode and the two endcap 
electrodes, d describes the dimension of the trap and is defined 


by d = 
electrode and 2zo the distance between the endcaps. The forces 
that the ion in the trap will suffer are: 


{i sci a 


A [222 + oe /2, where po is the inner radius of the ring 


where È, 5 -% - zê, and E, = fe - pêp. Therefore, the axial 
| motion is a harmonic oscillation parallel to the magnetic field 


with a frequency of 
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and the radial motion is a superposition of two eigenmotions 
with frequencies of 


We w? w? 
o= ZENT?” @) 
where 
ie LB. (4) 
m 


In order to remove unwanted ions and measure the ion’s 
mass, one method is to perform buffer gas cooling and rf ex- 
citations on the ion. Thus the equations of ion motion become 
more complicated than the ideal case mentioned above. 

A dipole field is created by an rf voltage with amplitude V4 
and frequency w,r applied with a phase difference of 180° be- 
tween two opposite segments of the ring electrode, and it gives 
an additional electric field [11, 12], for example, for the radial 
x-component: 


V, 
E, = costwst — rf): Ex, (5) 


where a is the inner radius of the trap. An azimuthal quadrupo- 
lar rf field with amplitude V} and frequency w,,¢ applied with 
180° phase shifts on sets of ring-electrode segments perpendic- 
ular to each other gives an additional electric field [10]: 
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where &, and Ê are the unit vectors on the x and y axes, respec- 
tively. 

Zhu et al. [13] studied the energy limitation for models to 
simulate the buffer gas cooling and suggested that the viscous 
drag force model should be used when the ion’s energy is less 
than ~5 eV/u. Thus the effect of the buffer gas on the ion motion 
in a Penning trap can be described as F = —6mv, where 6 is 
the damping parameter describing the effect of the buffer gas. 
With the ion mobility Ko, 6 can be written as 6 = TA ve : 
Here, g/m is the ion’s charge-to-mass ratio and P and T are the 
gas pressure and temperature in units of normal pressure and 
temperature Py and Ty [14], respectively. 

Thus the dynamical equations including all effects from the 
dipolar rf excitation and buffer gas cooling can be described as: 
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and those for the quadrupolar excitation and buffer gas cooling: 
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Here, the first part describes the ion motion in an ideal Penning 
trap, the second is from the dipolar/quadrupolar rf electric field 
and the third is from the buffer gas. 

As these are nonlinear differential equations, an analytic so- 
lution is difficult to obtain, and the method which has been 
known is using variable substitution twice VO = O) - 
wzBlt) X &, A (A) = VETO) and neglecting the high 
frequency terms [8, 10]. Here we obtain the numerical solution 
of the above equation by using Runge-Kutta method. 


3. Numerical solution by the Runge-Kutta method 


Common fourth-order Runge-Kutta method is used to obtain 
the numerical solution of ion motions in a Penning trap. Spec- 
ifying the form of first order differential equation as y’(t) = 
f(t,y) and defining its initial value as y(to) = yo, the follow- 
ing equations are introduced: 
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Table 1: Parameters used in the calculation 
h(s) B(T) Uo(V) d(mm) a(mm) m(u) P(Pa) 
1x10% 6.95 100 25.84 16 200 5 
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Figure 2: Axial ion trajectories in a plane along the magnetic field as a function 
of elapsed time. (a) without buffer gas; (b) with helium buffer gas at a pressure 
of 5 Pa. 


where A is the time interval, yg+ı is the approximation of y(tk+1) 
and 
fi = fte. yk) 
f= Fle + 3.ye + 3A (10) 
fa = fte + 5, Ye + 3f) 
fa = f(t + h, yk + hfs). 


Studies show that the error per step is on the order of h*, while 
the total accumulated error has order h* [15, 16]. Four slopes 
corresponding to each dependent variable need to be calculated 
by the same method mentioned above. The equations (7) and 
(8) can be rewritten as first order differential equations by defin- 
ing new variables x1 = x and yl = y. 

Table 1 shows some parameters used in the calculation. A 
proper value of interval A was chosen to reduce the accumulated 
error and running time for the calculation. The parameters B, 
U, d and a in the LPT [6, 17] are used. In order to show the 
effect of excitation clearly, larger values have been chosen for 
Va, the driving rf amplitude for dipolar excitation, V}, that for 
quadrupolar excitation and P, the helium buffer gas pressure 
for cooling. The trajectories of ion with mass 200 have been 
obtained for many different cases. 


3.1. Axial motion 


Figure 2(a) shows the ion trajectory in axial direction along 
the magnetic field without the presence of buffer gas. It is a har- 
monic oscillation with frequency w, and can be obtained from 
the analytic solution of equation (1). But with the presence 
of buffer gas, it is difficult to get the analytic solution with- 
out approximation. The numerical solution (figure 2(b)) shows 
that the ion will experience a damped vibration and the distance 
from the center of the trap will get smaller and smaller. Finally, 
the ions will be concentrated very close to the center of the trap. 


3.2. Radial motion without excitation 


In an ideal trap without the electric field for excitations and 
without the presence of buffer gas, the radial motion is a super- 
position of reduced cyclotron motion with frequency w and 
magnetron motion with frequency w_. The ion trajectory is 
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Figure 3: Radial ion trajectories in a plane perpendicular to the magnetic field 
without excitation. The trap center is at the point (0, 0) marked as a cross, and 
the ion starts at the point (1.0 mm, 0.0 mm). (a) without buffer gas; (b) with 
helium buffer gas at a pressure of 5 Pa. Note that the ion moves in a clockwise 
direction. 


shown in figure 3(a). The radius of cyclotron motion depends 
on the initial velocity of the ion, the greater the initial velocity, 
the larger the cyclotron radius. The radius of magnetron motion 
depends on the initial distance from the trap center, and a further 
distance makes a lager radius of the magnetron motion. Without 
the presence of buffer gas, the cyclotron and magnetron motions 
of the ion will keep its initial status. But with the buffer gas, the 
cyclotron motion will be damped and the magnetron motion 
will increase steadily but slowly, as shown in figure 3(b). 


3.3. Radial motion with azimuthal dipolar excitation 


The phase of the exciting dipolar rf field with respect to the 
initial magnetron motion has a very strong effect on the mag- 
netron radius [11, 12]. Depending on this phase difference one 
may actually decrease the magnetron radius. In all the calcula- 
tions shown in this section, we choose the phase difference as 
the most favorable value to enhance the magnetron radius. 

Figure 4(a) and (b) show the ion trajectories in radial di- 
rection when the azimuthal dipolar excitation is applied at 
w,f = w., for the cases without and with buffer gas cooling, 
respectively. The buffer gas has no obvious effect on the mag- 
netron motion and its amplitude increases continuously. More- 
over, all ions can be driven to a larger radius by this azimuthal 
dipolar excitation, because the first-order approximation of w- 
is independent of the ion’s mass. For the cyclotron motion, its 
amplitude becomes smaller and smaller when the buffer gas is 
present. 

Figure 4(c) and (d) show the similar results as figure 4(a) and 
(b), but with the applied dipolar excitation at w, = w+. Clearly, 
the dipolar field drives the cyclotron motion and the magnetron 
motion keeps constant without the presence of buffer gas. But 
with the buffer gas, the amplitude of the cyclotron motion is 
damped and that of the magnetron motion increases slowly and 
continuously. 

When the dipolar excitation is applied at frequencies other 
than w_ and w,, the picture is not very simple because other 
sidebands, for example, 2w_ and 2w,, may occur. However, if 
the excitation is performed at frequencies other than those par- 
ticular values, the field doesn’t drive the cyclotron or magnetron 
motion and it almost has no effect on the ion trajectories. The 
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Figure 4: Same as figure 3, but with azimuthal dipolar excitation. The left 
panels are obtained without buffer gas, and the right panels are with helium 
buffer gas at a pressure of 5 Pa. (a) and (b): w;f = w_; (c) and (d): wf = w+; 
(e) and (f): wf # w-, w, or other sidebands. Note that the ion moves in a 
clockwise direction. 


ion trajectories shown in figure 4(e) and (f) are very similar to 
those (figure 3) obtained without any excitation. 


3.4. Radial motion with azimuthal qradrupolar excitation 


Figure 5 shows the ion trajectories in radial direction when 
the azimuthal quadrupolar excitation is applied at the resonant 
mode, w,f = We = w, + w, and at the non-resonant mode, 
wrf # We, for the cases without and with buffer gas cooling. 
Initially there is only magnetron motion. 

Without the presence of buffer gas, the two eigenmotions are 
coupled and a full periodic conversion between the two motions 
is obtained after a certain time (depending on the amplitude V} 
of the excitation). For non-resonant excitation, the conversion 
is not complete. By scanning the frequency of this quadrupolar 
excitation, the coupled cyclotron frequency can be obtained by 
using the time-of-flight ion cyclotron resonance (TOF-ICR) de- 
tection technique [18] and then the mass of the ion can be calcu- 
lated by using equation (4). This is typically used for the mass 
measurement in most of the Penning traps around the world. 


When the buffer gas is present in the trap, the behaviour of 
the ion becomes quite different. When in resonance, if reason- 
able buffer gas pressure is applied, both the cyclotron and mag- 
netron motions decrease as a function of time and eventually 
the ion will go back to the center of the trap. On the contrary, 
when out of resonance, the ion will keep moving at the larger 
orbit and never go back to the center. Thus, this mass-selective 
centering can be used to purify the ion species and the con- 
taminations will be lost by using a very small pore during the 
transportation of the ions. 


In fact, when in resonance, the decrease of the motional am- 
plitudes as a function of time depends heavily on the damping 
constant. This has not been studied in detail and thus over- 
looked in the previous research papers. Figure 6 shows the ion 
trajectories in radial direction when the azimuthal quadrupolar 
excitation is applied at the resonant mode with different buffer 
gas pressures. When the buffer gas pressure in the trap is too 
low (Figure 6(a)), some periodic conversion could still occur as 
the gas can not cool the ion fast enough. As the pressure goes 
higher and higher, the ion can be cooled faster and faster, and 
both the cyclotron and magnetron motions then decrease as a 
function of time and eventually the ion will go back to the cen- 
ter of the trap (Figure 6(c)). But if the buffer gas pressure is too 
high then the increase in magnetron radius becomes significant, 
and the ion will never go back to the center (Figure 6(e)). So 
if one wants to use this method to recenter the ion of interest 
and purify the ion species, one has to set a reasonable buffer 
gas pressure. 


Figure 7 shows the centering time as a function of helium 
buffer gas pressure P and driving rf amplitude V, for azimuthal 
quadrupolar excitation. All excitations are performed at w,f = 
wc and all ions start at the point (1.0 mm, 0.0 mm). The cen- 
tering time is defined as the time-of-flight when the distance 
between the ion’s position and the trap center has decreased to 
less than 0.1 mm for 0.1 ms, which we consider the ion has 
been recentered. Obviously, to obtain a shorter centering time, 
one has to find a reasonable combination of P and V,. In our 
calculation, when different buffer gases, such as helium, neon, 
and argon, are used, the corresponding shortest centering times 
have been obtained to be 0.36 ms at ~10 Pa and ~28 V, 0.53 ms 
at ~5 Pa and ~28 V, and 0.90 ms at ~0.5 Pa and ~16 V, re- 
spectively. As heavier buffer gas is used, the optimal gas pres- 
sure becomes lower, maybe due to the increase of the stopping 
power. For a lower buffer gas pressure, the centering time be- 
comes longer, but too much buffer gas makes the ion move to 
the larger radius. For the driving rf field, its amplitude also 
affects the centering time, but we can not apply a very high 
amplitude practically in a real experiment, since it distorts the 
electrical field in the trap too much. So, to obtain a reasonable 
centering time in a real experiment, we recommend that the he- 
lium gas pressure should be set in the range of 0.01 — 5 Pa and 
the driving rf amplitude for quadrupolar excitation should be 
set in the range of 2 — 15 V. 
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Figure 5: Same as figure 3, but with azimuthal quadrupolar excitation. Initially there is only magnetron motion. The left and middle panels are obtained without 
buffer gas, and the right panels are with helium buffer gas at a pressure of 5 Pa. (a) and (b): w;f = We; (c) and (d): wf + We. (al) and (a2) show the first and second 
half of the conversion. Note that the ion moves in a clockwise direction. 
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Figure 6: Same as figure 3, but with azimuthal quadrupolar excitation at w,¢ = wc and with helium buffer gas at pressures of (a) 0.05 Pa, (b) 0.5 Pa, (c) 5 Pa, (d) 
30 Pa and (e) 50 Pa. Initially there is only magnetron motion. (al)—(a3) and (b1)-(b3) show successively the ion trajectories with helium buffer gas at pressures of 
0.05 Pa and 0.5 Pa, respectively. Note that the ion moves in a clockwise direction. 
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Figure 7: (color online) The centering time T as functions of helium buffer gas 
pressure P and driving rf amplitude V, for azimuthal quadrupolar excitation. 
All excitations are performed at w,f = We. 


4. Summary 


Penning traps have become very accurate tools for mass de- 
termination both on stable and unstable isotopes. In order to 
manipulate the ions in the trap, azimuthal dipolar and quadrupo- 
lar rf excitations and buffer gas cooling have to be applied, thus 
the ion motions become very complicated and the analytic solu- 
tions are impossible to obtain without approximation. We have 
studied in detail the numerical solutions of the ion motions by 
using the Runge-Kutta method for the different cases which can 
be used in the experiments for high precision mass measure- 
ments. The ion trajectories in the trap for those different cases 
have been obtained and summarized in graphical form. 

Ion motion in the trap depends heavily on the multipolarity of 
the rf field, its frequency and the buffer gas pressure. Azimuthal 
dipolar excitation at frequency w,f = w- can be used to drive 
all the ions to a larger radius. Azimuthal quadrupolar excitation 
at frequency w,f = we without the presence of buffer gas in the 
trap is generally used to measure the mass of the ion. For the 
recentering of the ion of interest and to perform the purification 
of the ion species, one has to set a reasonable buffer gas pressure 
in the trap and apply quadrupolar excitation at frequency w,f = 
We. 
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